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A Metal-Disordered Mott insulator-Metal heterostructure is studied at half-fiiling using unre- 
stricted Hartree Fock method. The corresponding clean system has been shown to be an insulator 
for any finite on site correlation. Interestingly we find that introduction of explicit disorder induces 
a metal-insulator transition at a critical value of disorder. The critical value corresponds to the 
point at which disorder nullifies the effect of onsite correlation. The wavefunctions are found to 
delocalize by increasing disorder, rendering the system metallic. 
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INTRODUCTION 

The field of hcttcrostructurcs, has seen tremendous 
interest in recent times both experimentally and 
theoretically Heterostructures fabricated out of 

strongly correlated systems, in particular have shown 
very novel properties. These emergent properties are 
sometimes very different from the properties of the in- 
dividual constituent materials. For example the het- 
erostructure of a Mott insulator and a band insulator has 
been shown to have metallic behaviour [lj. This is due to 
aggregation of delocalized electrons at the interface of 
such a heterostructure. On lowering the temperature, 
the system can even become superconducting [2[. 

In this paper, we investigate a metal-disordered Mott 
insulator heterostructure. The corresponding clean case 
system has recently received lot of theoretical attention 
and has been variously claimed to be a fragile Fermi 
liquid which shows metallic character at any finite cor- 
relation strength, and a " de facto impenetrable Mott in- 
sulator" @, with highly decaying spectral weight at the 
Fermi energy. In a recent work weQ have shown that 
the clean system is insulating for all finite values of cor- 
relation that we could work with. The system was also 
found to have multiple gaps in the energy spectrum, as 
a result of charge reconstruction on varying correlations. 

An interesting relevant question in this context is the 
role of disorder on the whole phenomenon. Experimen- 
tally heterostructures are fabricated with a certain degree 
of disorder present in them. The disorder can arise dueto 
non trivial local atomic and electronic structure due to 
the presence of dangling bonds and incomplete atomic 
coordination. With these fundamental issues in mind, 
we have investigated a heterostructure comprising of a 
metal-disordered Mott insulator-metal sandwich. Very 
interestingly, we find that the introduction of disorder 
can generate a metal insulator transition in the region 
where the corresponding clean system was gapped. This 



is because disorder is able to nullify the effect of correla- 
tion in a certain range of values and also delocalize the 
wave functions. The route to metallicity has very close 
connection to the metallic phase that is stabilised by in- 
creasing disorder from a clean Mott insulator, due to the 
closing of the Mott gap as shown by Aguiar et. al [§], 
using generalized DMFT methods. 

In the system we study here, there is charge reconstruc- 
tion between the planes and consequent inhomogenity 
in the underlying charge landscape. This inhomogenity 
was found to be sufficient to render the system insulating 
even in the gapless state of the clean system. Here we 
show that disorder is able to negate the effect of correla- 
tion over a small range of values and generates a metallic 
phase. 

We work with correlated disorder of the Thue-Morsc 
and Fibonacci type. Correlated disorder has been stud- 
ied extensively before U Ht ■ Theoretical studies of 
the quasi-periodic systems have shown that in the non- 
interacting limit the wave function shows a power law 
localization both in one[3l and twoflol] dimensions. A 
very recent experiment in a ID quasi-periodic optical 
lattice where the system was described by an Aubry- 
Andrc Hamiltonian, has shown exponentially localized 
states(Anderson localization) in the large disorder limit. 
A theoretical explanantion of the above result has been 
attempted using Fibonacci sequence in ID [l6[ 

Correlated disorder introduces a certain type of disor- 
der which docs not depend on the realisation, but still 
breaks translational sysmmetry. Its relevance in captur- 
ing the effect of interplay between correlation and disor- 
der has been shown recently by us[17|. The two chosen 
sequences have been wrapped over the quasi 2D system 
in the Mott layers. 

We use unrestricted Hartree Fock to study the problem 
as it is able to capture the role of disorder exactly and 
the effect of spatial variation in the order parameters 
both in plane and perpendicular to them. The effect of 
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correlation is captured in a fully self consistent manner. 
It works very well in higher dimensions at half filling 
and even in 2D [l8| and ID [l9[ gives results comparable 
to QMC and real space renormalization calculations. It 
cannot however capture the effect of large effective mass 
renormalization and the evolution of the sharp metallic 
peak in the spectral function at the Fermi energy. 



Model and Method 

We now present and discuss our work. We take the bar- 
rier of our metal-barricr-metal sandwich to be described 
by the single-orbital disordered Hubbard model, and the 
metallic planes by the tight-binding model. The Hamil- 
tonian for the system is 

ft = ~ ti Uia ~ X! 

i<j ijaa 

+ U «( T Ho l f - h(n la i - r)- (1) 



Here i,j is the site index. The operator cj acr (c icr ) cre- 
ates (destroys) an electron of spin a at site i. We set both 
the in-plane and the inter plane hopping to be nearest 
neighbor only, and equal to t. The lattice structure is 
that of a simple cubic lattice, q denotes the site poten- 
tial at the ith site which follows either a Thue-Morse or 
a Fibonacci sequence in the barrier planes and are zero 
in the metallic planes. 

A Thue-Morse sequence is generated by following the 
rule: A — > AB and B — > BA, where A and B are the 
two different sites. In our case site B has a higher site 
potential than site A {ts > £a, W = 6b — £a, being the 
disorder strength) . Thue-Morse chain in a particular gen- 
eration looks like: ABB AB AABB AAB ABB A.. . The A 
and B sites are found in equal numbers in any Thue 
Morse sequence of even entries. The generating scheme 
of Fibonacci sequence is A — > AB and B — > A. A typi- 
cal Fibonacci chain in a particular generation looks like: 

ABA ABABA ABA ABABA ABABA Here sites A 

and B follows a golden mean ratio. These two sequences 
correspond to two different impurity site concentration in 
the system. In Thue Morse 50% of the barrier sites are 
B type, while in Fibonacci the corresponding number is 
38% . It is therefore expected that the effect of disorder 
would be more strong for the Thue Morse disorder com- 
pared to the Fibonacci disorder. These sequences have 
been wrapped over the quasi 2D lattice in the barrier 
planes. The interaction term is written in a particle-hole 
symmetric fashion. The chemical potential [i is calcu- 
lated by demanding that there be exactly N electrons in 
the problem and is calculated by taking the average of 
the N/2 th and the (N/2 + 1) th energy level. We label 
the planes by a with a = m or a = 1 corresponding to 
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FIG. 1: a. Variation of gap for [7 = 5 with increasing disor- 
der for Thue Morse and Fibonacci sequences, b. IPR with 
increasing disorder for U — 5 for Fibonacci and Thue Morse 
sequences 



the metallic layers. All the other a values in between 
correspond to the barrier planes. Thus U\ = U rn = 0, 
while for all other a, U a = U. 

In a recent work 3 we had studied the interface of 
a quasi 2D barrier with onsite correlations sandwiched 
between two metallic planes and shown that even a very 
small correlation stength renders the system insulating. 
The system for very low U behaves like a gaplcss insu- 
lator, while for larger U, the system opens a gap at half 
filling and then behaves like a gapped insulator. For even 
higher U, the system shows the presence of multiple gaps. 
In this work we will fix U to be in the region where the 
system has opened a single gap. We then turn on the 
disorder strength from zero and study the evolution of 
this sytem. 

We use unrestricted Hartree Fock method to address 
the problem as we do not want to restrict or bias the 
solutions in any particular manner. The price we pay in 
the process is the time taken to attain self consistency 
in the solutions, upto the high precision that we seek. 
This precision value is 10 -7 in the square of the den- 
sity of both up and down electron density. We need to 
perform about 100 iterations away from the transition 
point to reach self consistency, while close to the metal- 
insulator transition point we need to perform more than 
1000 iterations to reach the same consistency. For disor- 
der strength even higher than the transition value, it is 
very difficult to reach self consistency as the numerical 
solution vacillates between several close lying minima so- 
lutions. Hence we avoid this region of parameters in this 
work. The initial seed that we provide to the system be- 
comes increasingly important as the system approaches 
the metal insulator transition point. Deep in the Mott 
insulator regime we find that though the system is not 
highly sensitive to the initial seed in determining the out- 
come of the final ground state, a mixed antiferromagnetic 
seed gives the lowest energy for both Fibonacci and Thue 
Morse type of disorder. Close to the transition point we 
find the initial seed that gives the lowest energy is differ- 
ent for the two types of disorder. This is because close 
to the transition point the effective single particle poten- 
tial landscape that emerges due to competition between 
correlation and disorder becomes highly sensitive to the 
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FIG. 2: Spin and charge profiles for Fibonacci type of disor- 
der, a. Spin at each site for U — 5 and W = 0.6. b. Charge 
at each site for U — 5 and W = 0.6. c. Spin at each site for 
(7 = 5 and W = 0.82. d. Charge at each site for U — 5 and 
W = 0.82. 



concentration and distribution of the impurity sites. 



Calculated quantities 

We study the effect of increasing W on the gap at half 
filling, the spin and charge profile and the Kubo conduc- 
tivity at zero temperature. The charge at a particular 
site is simply calculated as C = + ?&i 4.. The spin at 
a particular site is given by S = (rii.f — "1,4.) |- The op- 
tical conductivity is calculated using the Kubo formula, 
which at any temperature is given by: 



a, p 



with A = ire /hao, ao being the lattice spacing, and 
n a = Fermi function with energy e a — fi. The f a p 
are matrix elements of the current operator j x = 
itJ2it <j( c l+xa <7 C i,v ™ h.c), between exact single particle 
cigenstatcs \ip a ), \ipp) t etc, and e a , ep are the correspond- 
ing eigenvalues. In this paper, conductivity /conductance 
is expressed in units of A = ire 2 /hao. 

We calculate the 'average' conductivity over 3-4 small 
frequency intervals Aoj(Aw = nuj r , n = 1,2,3,4), and then 
differentiate the integrated conductivity to get <t{u) at 
lu = nui r , n = 1.2.3.4 [rH, l20j |. The uj r is taken to be twice 
the lowest frequency that can be accessed of a finite size 
lattice (uj r = 2 * D/N), where D is the bandwidth and N 
is the total number of sites in the system. 

For the quasi 2D geometry taken by us N = m x L 2 , 
where m is the total number of layers (metal + insulator) 
and L is the total number of sites along both x and y 
direction, in the numerical work. 



Results and Analysis 

Our results in this paper corresponds to U — 5 in pres- 
ence of disorder. We have chosen this value as represen- 
tative of the region where the clean system is a Mott 
insulator. As we increase the disorder strength, the gap 
at half filling closes and the system becomes metallic. 
There is also a gradual derealization of the wavefunction 
as the disorder strength is increased. The spin order pa- 
rameter becomes highly dispersed in each of the planes as 
the metallic phase is approached. The peak conductivity 
values are roughly the same for both types of correlated 
disorder considered in this work. 

Fig la shows the phenomenon of the gap at half fill- 
ing, due to correlation, closing due to the presence of 
disorder. Disorder introduces states at half filling. The 
concentration of B sites is more for the Thue Morse se- 
quence compared to the Fibonacci sequence. Hence the 
Thue Morse sequence is able to kill the gap more quickly 
than the Fibonacci sequence. Thus we see that while the 
Thue Morse sequence kills the gap to a value less than 
the finite size spacing at roughly W = 0.6, the Fibonacci 
sequence is able to achieve the same reduction in gap at 
around W = 0.82. 

Fig lb looks at the change in the extent of localization 
of the wavefunction as we increase the disorder through 
the Inverse Participation Ratio(IPR). High IPR values 
indicate higher localization while low IPR values indicate 
delocalized wavefunctions. The IPR decreases monoton- 
ically with increasing disorder, thus indicating that the 
wave functions become more delocalized. The calcula- 
tion has been done for both Thue Morse and Fibonacci 
sequence of disorder. The derealization process is more 
rapid for the Thue Morse type and is more gradual for 
the Fibonacci sequence. We have not been able to go 
beyond the transition point due to the onset of a glassy 
phase which makes it very difficult to calculate the cor- 
rect configuration accurately. 

From our previous studies Q we know that for m = 5 
which means 3 intervening Mott layers, there is a gap 
at half filling for U = 5 accompanied by the induction 
of anti ferromagnetic ordering. The Mott layers induces 
an anti ferromagnetic ordering even in the metal planes. 
The ordering is non uniform as one moves along the z-axis 
with the maximum value at the central plane. 

Fig 2 and Fig 3 shows the spin and total charge at 
each site for Fibonacci and Thue Morse type of disorder 
respectively, for L = 22, m = 5, U = 5. Thus in the 
figures, the first L 2 sites correspond to the lowest metallic 
plane, the next L 2 sites correspond to the adjacent Mott 
plane, the next 1? to the central layer Mott plane and so 
on. There is obvious symmetry about the central plane 
in the spin and charge order parameters. 

Fig 2a and 2c show the spin profile for the 
insulating(VF = 0.6) and metallic (W = 0.82) regime re- 
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CONCLUSION 



A metal-disordered Mott insulator-metal interface het- 
erostructure is studied. A disorder induced metal insula- 
tor transition is observed, over a small range of disorder 
values, in the range where the clean system is gapped. 

Sanjay Gupta would like to thank the DST for provid- 
ing the grant under the project "Electronic states and 
transport in mesoscopic/nanoscopic systems". 



FIG. 3: Spin and charge profiles for Thue Morse disorder, a. 
Spin at each site for {7 = 5 and W = 0.5. b. Charge at each 
site for U — 5 and W = 0.5 c. Spin at each site for U — 5 and 
W = 0.63, d. Charge at each site for U = 5 and W = 0.63 
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FIG. 4: a. Insulator-Metal-Insulator transition observed with 
increasing disorder strength of the correlated Fibonacci type. 



spectively, while Figs 2b and 2d show the charge profile 
for the same values of W. Figs 3a and 3c shows the cor- 
responding spin profile for the insulating (W — 0.5) and 
metallic (W = 0.63) regime respectively, for the Thue 
Morse sequence, while Figs 3b and 3d show the charge 
profile for the same values of W. 

We find disorder increases the total charge buildup in 
the metallic layers. The charge depletion from the Mott 
layers just adjacent to the metallic planes and the the 
charge accumulation in the metallic planes increases with 
increasing W . The spin profile becomes totally dispersed 
for the metallic phase in contrast to a less dispersed pro- 
file in the insulating phase. 

From the Kubo conductivity plots shown in Figs 4a 
and 4b wc find that as we increase disorder for fixed U = 
5, the system shows a metal insulator transition(MIT). 
The conductivity value changes non monotonically with 
increasing disorder in the metallic phase. 



t 

[1 

[2 
[3; 
[4 
[5 

[6 

It 

[8 
[9 

[10 

[11 

[12 



[13 

[14 
[15 

[16 
[17 

[18 
[19 
[20 



sanjayl@bose.res7in] 



trbkrm01@yahoo.co.in 



A. Ohtomo et al., Nature (London), 419, 378 (2002), 
Nature (London), 427, 423 (2004). 

S. Thiel, G. Hammerl, A. Schmehl, C. W. Schneider, and 

J. Mannhart, Science, 313, 1942 (2006).. 

Kenji Yonemitsu, Nobuya Maeshima and Tatsuo 

Hasegawa, Phys. Rev. £.,76, 235118, (2007) 

W. C. Lee and A. H. MacDonald, Phys. Rev. £,74, 

075106 (2006). 

S. Okamoto and A. J. Millis, Nature (London),428, 630 
(2004); Phys. Rev. B, 70, 241104(R) (2004); Phys. Rev. 
£,72, 235108 (2005). 

R. W. Helmes, T. A. Costi, and A. Rosch, Phys. Rev. 
Lett. 101, 066802 (2008). 

H. Zenia, J. K. Freericks, H. R. Krishnamurthy, Th. Pr- 
uschke , Phys. Rev. Lett, 103, 116402, (2009). 
Sanjay Gupta, Tribikram Gupta arXiv 1003.4642 
Aguiar M. C. O. et. al. , Phys. Rev. Lett. 102, 156402 
(2009) 

Merlin R. et. sl.,Phys. Rev. Lett, 55, (1985) 1768. 
Delahaye J. et. &l.,Phys. Rev. B, 67 (2003) 214201. 
Sanjay Gupta et. al., Physica B, 355 (2005) 299; De Brito 
P. E. et. al. Phys. Rev. B, 73 (2006) 014301; Mauriz P. 
W. et. al., Phys. Rev. B, 63 (2001) 184203; Vasconcelos 
M. S. et. al., Phys. Rev. B, 59 (1999) 11128. 
Roati G. et. al. , Nature, 453 (2008) 895. 
Kohomoto Mahito et. al. ,Phys. Rev. £,34 (1986) 563; 
Sutherland B.,Phys. Rev. B, 34 (1986) 3904, Phys. Rev. 

B, 39 (1989)5834. 

Modugno M., New J. Phys., 11 (2009) 033023. 
Tribikram Gupta and Sanjay Gupta, Euro Physics Let- 
ters 88, 17006 (2009). 

Lin H. Q. and Hirsch J. E., Phys. Rev. B 35, 3359 (1987) 
Sanjay Gupta et. al. Phys. Rev. B 63, 125113 (2001) 
Sanjeev Kumar and Pinaki Majumdar, Europhys. Lett 
65, 75 (2004) 



